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Abstract. - We present an inelastic neutron scattering investigation of Li2Cu02 detecting the 
long sought quasi-lD magnetic excitations with a large dispersion along the Cu02-chains studied 
up to 25 meV. The total dispersion is governed by a surprisingly large ferromagnetic (FM) nearest- 
neighbor exchange integral Ji = —228 K. An anomalous quartic dispersion near the zone center 
and a pronounced minimum near (0,0.11,0.5) r.l.u. (corresponding to a spiral excitation with 
a pitch angle about 41°) point to the vicinity of a 3D FM-spiral critical point. The leading 
exchange couplings are obtained applying standard linear spin-wave theory. The 2"*^ neighbor 
inter-chain interaction suppresses a spiral state and drives the FM in-chain ordering below the 
Neel temperature. The obtained exchange parameters are in agreement with the results for a 
realistic five-band extended Hubbard Cu 3d O 2p model and LSDA-|-(7 predictions. 



. 1. INTRODUCTION. - Li2Cu02 is the first [1] 
dnd the most frequently studied compound of the growing 
class of edge-shared spin-chain cuprates [2-7]. Owing to 
its structural simplicity with ideally planar Cu02 chains 
(isee Fig. 1) it has been considered as a model quasi-one- 
dimensional (ID) frustrated quantum spin system. In al- 
most all edge-shared cuprate chain compounds the spins 
are expected to be coupled along the chains via near- 
est neighbor (NN) ferromagnetic (FM) and next-nearest 
neighbor (NNN) antiferromagnetic (AFM) exchange in- 
teractions, Ji and J2, respectively. Due to the induced 
frustration FM and spiral in-chain correlations are com- 
peting. While in the ID model the ground state is gov- 
erned by the ratio a = — J2/J1, the actual 3D magnetic 
order sensitively depends on the strength of the inter-chain 
couplings and anisotropy. In Li2Cu02, below « 9K 
[8, 9] a long-range coUinear commensurate AFM inter- 
chain with FM in-chain (CC-AFM-FM) magnetic ordering 



evolves. However, a proper understanding, necessary for 
a critical evaluation of theoretical studies, especially of 
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Figure 1: (Color online) Left: The crystallographic structure 
of Li2Cu02 comprises two AFM coupled CUO2 spin-chains per 
unit cell running along the 6-axis (orange • - Cu^+, red • 
- O^^, bright blue - Li"^). The unit cell is indicated by 
the outer black cuboid. Right: the main intra- and inter-chain 
exchange paths, Ji , J2 , and J2 marked by blue arcs and dashed 
lines, respectively. Notice the frustration introduced by an 
AFM inter-chain coupling for any non-FM in-chain ordering. 
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Figure 2: (Color online) (a) Constant energy scans for momentum transfer along the chains with q = (0 if 1) at IN12 
(Ai5 < 3meV; open symbols) and q = (0 1 + A" 0) for INS data (Ai5 > 5meV; filled symbols), (b) Intensity map of the 
low energy spectrum, interpolated from the IN12 data points (•). The sizeable intensity below the gap energy, i.e. at energy 
transfer < 1 meV and qh around -0.05 is a spurious Bragg tail, (c) The dispersion perpendicular to the chains was measured 
by constant Q-scans; lines are GAUSSian fits. Filled symbols: the measured intensity at the zone center for T > Tn. 



electronic/magnetic structure calculations [10-14], is still 
missing. In particular, this concerns a precisely enough 
knowledge of the main exchange interactions. The knowl- 
edge of realistic values is helpful also for the understand- 
ing of related "frustrated ferromagnets" [15, 16] such as 
CaaYaCugOio [17], La6Ca8Cu2404i [2] with FM in-chain 
ordering and LiVCu04, LiCu202 with helimagnetism and 
multiferroicity, all being of considerable current interest. 

A previous inelastic neutron scattering (INS) study 
aimed to determine the exchange integrals was not con- 
clusive [18]. It revealed an anomalous low- lying branch 
of hardly dispersive and overdamped spin excitations in 
chain direction. Its linear spin-wave (LSW) analysis re- 
sults in unrealistically small in-chain exchange integrals. 
The missing, but expected, dispersive quasi-lD spin chain 
excitation remained as a challenging puzzle for the com- 
munity [10,11,19]. In Sec. 3 we present new INS data 
which unambiguously show the presence of a strongly dis- 
persive in-chain spin mode. The main exchange integrals 
are derived applying the LSW-theory [20] and in Sec. 4 we 
compare our results with those of related chain cuprates as 



well as with predictions of band structure and cluster cal- 
culations. A criticism of improper CURIE- Weiss analysis 
of spin susceptibility data is provided and consequences 
for the direct FM Cu-0 exchange parameter Kpd entering 
extended Hubbard models are discussed. 

2. EXPERIMENTAL. - A single crystal of 
^Li2Cu02 was grown for INS experiments by the trav- 
elling solvent floating zone technique under high pres- 
sure [21]. In order to avoid vaporization of Li20 dur- 
ing growth a 4:1 Ar:02 atmosphere at 50 bar was chosen. 
A fast growth rate of lOmm/h inhibits growth of impu- 
rity phases. Isotope enriched ^Li was employed to avoid 
the significant neutron absorption coefficient of ^Li. The 
sample was characterized by X-ray powder diffraction, po- 
larized light microscopy, magnetization and specific heat 
measurements. By X-ray powder diffraction no impurity 
phase was found. The macroscopic magnetization and spe- 
cific heat data of the sample agree with literature data, 
i.e. AFM order is found below Tn = 9.2 K and a weak FM 
component evolves below T2 « 3K [9,22,23]. 
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INS experiments were performed with thermal and cold 
neutrons at the three-axis-spectrometers INS and IN12 at 
the ILL, Grenoble. France. Four single crystals with a 
total mass of 3.8 g were mounted together in the (0 K 
L) scattering plane with a resulting sample mosaicity of 
3°. For both instruments, focusing PG(002) monochro- 
mator and analyzer have been utilized. The measure- 
ments at IN8 were taken with fixed final momentum 
fcf = 2.662 A with PG-filtcr on kf. IN12 was configured 
with kf = 1.5 A and Be- filter on fcf. Most scans have 
been done in the CC-AFM-FM phase at T = 4.1 K well 
above another not yet well understood magnetic phase be- 
low T2 ■ Anyhow, the observed changes of the INS spectra 
in this low-T phase (not shown here) arc weak. 

3. RESULTS AND DATA ANALYSIS. The 

results of our INS studies are summarized in Fig. 2. Rep- 
resentative spectra of constant energy scans as taken at 
IN12 and IN8 for moment transfer along the chains (&*) 
are displayed in Fig. 2(a). The main result is the obser- 
vation of a highly dispersive excitation which is strong 
at the magnetic zone center and significantly weakens at 
higher energies. With the chosen experimental setup the 
magnetic branch could be traced up to energy transfers of 
25meV. The measured data points along (0 K 1) taken 
at INI 2 with cold neutrons are summarized in the color 
map Fig. 2(b). Note, that the reflections are periodic with 
the magnetic unit cell and their strongly reduced inten- 
sity above Tn observed up to energy transfers of 15mcV 
does confirm their magnetic nature (see Fig. 2 (c)). At 
the magnetic zone center (0 1), a gap of A = 1.36 meV 
is observed. The excitations for momentum transfer along 
(0 1 -|- L) are only weakly dispersive in agreement with 
the results of Ref. [18]. Respective constant q-scans are 
shown in Fig. 2(c). Note, that the mosaicity of the sample 
broadens the excitations along L which is less pronounced 
for moment transfer along K due to the longer 6*-axis. 

As shown in Fig. 2(b) wc observe further inelastic fea- 
tures for moment transfer along the chain. In addition, 
the data in Fig. 2(b) also exhibit weak and presumably in- 
commensurate (IC) magnetic scattering below the magnon 
gap energy. The origin of these low-energy excitations 
is not yet clear and will be addressed in future studies. 
Furthermore, there is a continuous feature appearing at 
double the energy of the anisotropy gap which is possibly 
attributed to two-magnon scattering. We also mention 
that we have observed low-lying and strongly broadened 
excitations along h* similarly as in Ref. [18]. However, the 
intensity of these excitations was roughly two orders of 
magnitude weaker than that reported ibidem. 

According to ESR measurements [24] the exchange in- 
teractions show an uniaxial anisotropy with the easy-axis 
directed along the crystallographic a-axis. We describe the 
corresponding Cu momenta by the spin-HAMiLTONian 

m,r 



where m enumerates the sites in the magnetic (Cu) lattice, 
the vector r connects sites with an exchange coupling Jr 
[25]. The z-axis is taken along the easy-axis, i.e. the a-axis. 
Within the LSW-theory [20] , the dispersion-law reads 

= \l {J^' - Jo' + Jo" - d)' - {J^^y , (2) 

where Jq = (1/2) Jr cxp («qr) is the Fourier transform 
of the in-chain exchange integrals, and analogously for the 
inter-chain integrals Jq. The exchange anisotropy D = 
Jq — Jq ^ ^ Jq ~^ Jq' causes the abovemcntioned spin gap 
A — llIq in our case (see Fig. 3). Their relation reads 



A^^D[D-2J^y), or D ^ - ^ [j^^) + . 

(3) 

In the summations over r we retain only the leading 
terms (see Fig. 1). According to LSDA-|-f7 based mag- 
netic structure calculations they are given by the fol- 
lowing in-chain integrals: Ji, J2, J3, (corresponding to 
r = b, 2b, 3b respectively) and inter-chain integrals: 
Jill, Ji3i, (corresponding to rm = (a + b + c) /2, riai = 
(a -I- 3b -I- c) /2). For or QcC = tt the inter-chain dis- 
persion caused by Jq^ vanishes and Eqn. (2) simplifies: 

= - Jo" + Jo' J^y + Al , (4) 

i.e. the single-chain dispersion can be read off directly by 
subtracting the effective gap Ai, only. 

Our INS data are well fitted by GAUSSian distribu- 
tions. The maxima of the main branch of the spec- 
trum were analyzed within LSW-theory Eqn. (2). The 
inspection of Fig. 3 reveals a strong in-chain dispersion 
and a much weaker one in the perpendicular c and a 
(not shown here) directions. The results of the fit are 
given in Tablet. The total width W = uj^r/b — A = 2 | 

Ji + J3 + ./5 + 0.5 (j^^ - i?) I -A « 2 I Ji I, i.e. 

the in-chain dispersion yields a direct measure of the NN 
coupling since W is wnaffected by J2, J4, Jg • • • . Supposing 
a monotonous behaviour of LUq up to the zone boundary 
at Q = 0.5 (corresponding to q = 7r/6), from the mea- 
sured part of the full Wg-curve one obtains already a rig- 
orous lower bound for | Ji [^ 150 K. Our analysis shows 
that the in-chain NN interaction Ji is strongly FM, but 
frustrated by an AFM NNN-coupling J2 which affects the 
shape of tOq. In comparison, the AFM inter-chain coupling 
is weak, clearly demonstrating the magnetically quasi-lD 
character of the compound. For the in-chain coupling we 
find a = — J2/J1 — 0.33, unambigously above the critical 
ratio ttcrit — 1/4 for an isotropic HEiSENBERG-chain [26]. 
Note that the dispersion near the zone center behaves like 
uj{q) oc to be discussed below. We confirm also theo- 
retical predictions [12, 14] (see also Tab. 1 for our results) 
that the main inter-chain coupling is indeed the NNN cou- 
pling along (|-, ^, |). The NN inter-chain exchange has 
a negligible effect on the dispersion in the (0 K i)-plane. 
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Figure 3: (Color online) (a)-Constant energy scans for the data points shown in part (b) (right panel) near the minimum. (b,c) 
- INS data points and LSW fits (red line) of the magnon dispersion (b): along q =(0 K 1), (0 1 + L), (0 K 1.5) for _E-transfer 
below 5meV, compared with Ref. [18] (black dotted line); (c) - along Q = (0 if 0) up to high i5-transfer. Solid red line: our fit 
(see Tab. 1). For comparison the LSW dispersions predicted in Ref. [19] (green dashed-dotted and dashed double dotted lines), 
by our L(S)DA+C/ (blue dashed line) calculation and the five-band Hubbard model (blue dotted line), both with the added 
experimental spin gap are shown, too. 



Although further couplings can not be accessed from our 
fits, the main J values given here do not change much, if 
the INS data are analyzed in more complex models with 
additional exchange paths, especially J3 and Ji. When 
taken into account, both remain small ('^ 4 K and 1 K, 
respectively) in full accord with L(S)DA+C/ (sec Tab. 1). 

4. DISCUSSION. - The actual CC-AFM-FM or- 
dering seemingly contradicts the IC "spiral" phase ex- 
pected for a frustration ratio of a « 0.33 in a ID-approach 
(in the sense of the wave vector go 7^ 0, 7r/6 where the mag- 
netic structure factor S{q) becomes maximal [26]). Hence, 
the obtained relatively small but frustrated AFM inter- 
chain coupling may hinder the spiral formation. Thus, 
the 3D critical point compared with is upshift ed: 

^3D,iso ^ ^ID (1 + + + 25/33 + •••), (5) 

where /3„ = — J„/Ji. Eqn. (5) has been derived in the 
isotropic (iso) case [27]. Ignoring all other very weak inter- 
chain couplings J„ we arrive with our results J2 = 9.04 K 
and Ji = -228 K at a^^^'™ = 0.339. Anisotropics (aniso), 
as found here, further stabihze the CC-AFM-FM state. 
From Eqns. (2,3) we estimate finally a^^^^"'"" w 0.39. The 
combined effect of AFM inter-chain coupling and easy- 
axis anisotropy is also responsible for the anomalous q^- 
dependence of the spin excitations mentioned above. In 
the limit qi,b ^ we expand Eqn. (2) and obtain 

Lo{q)^^ + Ar{qbhf+Br{qbh)\ (6) 



and the quadratic dispersion vanish exactly at 

«o = 7 + ^) = 0.33098, <5 = 1--^. (7) 

Accidentally Li2Cu02 is very close to this point and its 
dispersion near the zone center is quasi- quartic. But in 
the presence of a spin gap A caused by the anisotropy 
D this vanishing of the quadratic dispersion doesn't yet 
signal an instabihty of the CC-AFM-FM state. 

Near the Z-point (0, 0, tt/c) the quadratic coefficient Az 
is already essentially negative (see Fig. 3 (b, right panel)). 
Since along the line Z — R{0, n/b, tt/c) the inter-chain dis- 
persion vanishes one can easily read off the ID Fourier 
components of the exchange interactions Jq^ (see Eqn. 
(4)). A similar rare situation occurs in the 2D frus- 
trated CsCuCU system in a high magnetic field above 
its saturation limit [28]. The clearly visible minima at 
qtfib cos^i (l/4a) « ±0.72 = ±0.11 (r.l.u.) correspond 
to the two equivalent propagation vectors of a low-lying 
spiral excitation [29] with a pitch angle of about 41.2° 
above a CC-AFM-FM ground state observed to the best 
of our knowledge for the first time. 

Next, we briefly compare our results with those obtained 
so far by INS-studies for CazYzCusOio (CYCO) [17] with 
a similar FM in-chain ordering and a frustrating AFM 
inter-chain interaction. There the reported Ji read -80 K 
and -93 K for fits where J2 = and J2 = 4.6 K, re- 
spectively. However, such tiny values of J2 are unlikely 
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Table 1: The fitted exchange integrals (in K) as determined from the INS data using Eqns. (2,3) compared with microscopic 
theory (see text) and other recent theoretical results. Values in parentheses are estimates from less accurate fits. 





Ji 
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J2 = ~0L ■ Jl 


J3 


J2 


D 


Jl 


INS / present work 


-228 ±5 


0.332 ± 0.005 


76 ±2 


(3.8) 


9.04 ±0.05 


-3.29 ±0.2 


(1) 


3d02p 1 present work [37] 




-218 


0.30 


66 


-0.4 








3d02p [30] 




-143 


0.23 


33 


-1 








3(i02p [10] 




-103 


0.47 


49 


-2 








two-chain plicnonienol. [19] 




-100 


0.40 


40 




16 




16 


LSDA+L/, L/ = 6 eV [38] 




-216 ±2 


0.31 


66 ±2 


5±2 


13 ±2 




0±2 


GGA L/ 6 eV [14] 




-171 


0.60 


98 




18 




0.23 


RFPLO, LAPW+SO [40] 














-15.6 





[10,19,30]. For a standard Cu-0 hybridization a much 
larger value is expected [31] in accord with the observed 
sizable part of the total magnetic moment (22 % ) re- 
siding at O [32]. Re-fitting their INS data yields Ji,J2 
values of the same order as we found for Li2Cu02 (LCO) 



CYCO 



JLCO and jpco _ Q^^jhCOy ^ detailed com- 



parison of both systems will be given elsewhere [31]. With 
respect to their large Ji values the question may arise why 
they have been not recognised so far in analyzing thermo- 
dynamic properties? In this context a critical evaluation 
of the reported AFM " CuRiE- Weiss" (CW) temperatures 
©cw"^ ~ -15 K [33] or small FM values: 5-10 K [34] and 
©cwP ~ -40 [8, 18] or -8 K [35] is very instructive. All 
these data have been derived from the linear fits of inverse 
susceptibility plots below 300-400 K (of the type denoted 
as "pscudo"-CW-lincs in Fig. 5). But here we estimate 



e 



CW 



Jl I -J2 



^i7it — ch 



Jl + J 2 



> +54 K, 



where the inter-chain coordination number Zint-ch = 4, 8 
for CYCO and LCO, respectively. Thus, we arrive at 
about 60 K and 58±4 K, respectively. The inspection 
of Fig. 4 clearly shows that very high T [36], far above 
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Figure 4: (Color online) The inverse spin susceptibility as mea- 
sured (o) vs high T-series expansion in 1"* (correct CW-law) 
and 2"'' order. Dashed and dotted lines: "pseudo"-CW-laws. 



any available data would be required to extract 0cw 
from l/x(T)-data. To satisfy the incorrect 6cw-values 
strongly underestimated Ji -values have been adopted in 
Refs. [10,18,33-35]. Even the 2"^^ order of high T-series 
expansion (HTS) approaches the experimental curve only 
above 400 K. Hence, even more any attempt to detect even 
a CuRiE-law [10] near 300 K must fail. Our Ji-values 
for LCO and CYCO from INS-data provide support for 
the large value we found in Li2ZrCu04 from thermody- 
namic properties [6] but puzzle the tiny value reported for 
LiVCu04 [5] with a similar Cu-O-Cu bond angle. 

Finally, we turn to a microscopic analysis. In Tab. 1 and 
Fig. 3 we compare our INS derived exchange integrals with 
theoretical results. First, we list the in-chain couplings J„ 
as obtained from the mapping of a five-band extended 
Hubbard pd model (on open chain Cu„02ri-i-2-clusters 
n = 5, 6) [37] on a corresponding Ji- J2- J3-HEISENBERG 
model (see Fig. 1). But here, to reproduce the main ex- 
perimental exchange integrals, a refinement has been per- 
formed most importantly by considering a larger direct 
FM Cu-0 exchange Kpd = 81 mcV compared with 50 meV 
adopted in Ref. [19]. We note that practically only Ji 
is significantly affected by Kpd- Thereby | Ji |oc Kpd 
holds approximately. Notice that the contribution of Kpd 
is much more important for the large negative (FM) value 
of Jl than that of the intra-atomic FM Hund's rule cou- 
pling on O. Since the available spectroscopic data at 300 K 
depend only weakly on Kpd not much is known on its mag- 
nitude. In the past Kpd has been used mostly as a fitting 
parameter for thermodynamic properties ranging from 50 
to 110 mcV for CuGeOs [10,39]. The INS data reported 
here provide a unique way to restrict its value phcnomeno- 
logically and opens a door for systematic studies of this 
very important interaction and well-founded comparisons 
with other edge-shared Cu02 chain compounds. 

Secondly, in the LSDA+C/ there is practically only one 
adjustable parameter Ud — Jh, where Ud denotes the 
Coulomb onsite repulsion (between 6 and 10 eV) and Jh 
denotes the intra-atomic exchange (« 1 eV) both on Cu- 
sites. Comparing the total energy of various ordered mag- 
netic states, a set of in-chain and inter-chain integrals can 
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be derived [14]. As a result one arrives again at very close 
numbers to our INS derived set [38]. Noteworthy, both 
ED a well as the LSDA+C/ provides a justification to ne- 
glect any long-range exchange beyond the third NN. The 
latter also explains why there is only one important inter- 
chain exchange integral J2 {132) ■ The excellent agreement 
between the INS-data analyzed in the simple LSW-theory 
and the theoretical results / predictions suggests that in the 
present case the effect of quantum fluctuations as well as 
of spin-phonon interaction seems to be rather weak. The 
former point is also supported by the relatively large value 
of the magnetic moment m « 0.96/^b [8,9] in the ordered 
state below Tfq and a consequence of the fact that the FM 
state is an eigenstate of the ID spin-model in contrast to 
the Neel state. Concerning the value of the anisotropy, 
there is no good agreement between contemporary DFT 
calculations [40] and much smaller values obtained in var- 
ious experiments [18,24] including our data (see Tab. 1). 

5. SUMMARY. - The main results of our INS study 
are (i) the relatively large dispersion of spin excitations in 
the Cu02 chains of Li2Cu02 due to the large value of the 
FM NN in-chain coupling Ji and (ii) the observation of a 
low-energy spiral excitation over a commensurate collinear 
Neel ground state in the vicinity of the 3D critical point 
above the corresponding ID point. The obtained main ex- 
change integrals can be approximately reproduced adopt- 
ing an enhanced value for the direct FM exchange Kpd 
between Cu 3d and O 2p states within an extended five- 
band HuBBARD-model. Further support for the empirical 
exchange integrals comes from L(S)DA+C/ calculations, if 
a moderate value of U somewhat smaller than the Ud in 
exact diagonalization for the extended HuBBARD-model is 
employed. The achieved detailed knowledge of the main 
exchange couplings derived from the INS-data provides a 
good starting point for an improved general theoretical 
description of other Cu02-chain systems and to adress a 
microscopic theory of their exchange anisotropy. 
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